High-lying single-particle modes, chaos, correlational entropy, and doubling phase 

transition 
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Highly-excited single-particle states in nuclei are coupled with the excitations of a more com- 
plex character, first of all with collective phonon-like modes of the core. In the framework of the 
quasiparticle-phonon model we consider the structure of resulting complex configurations using the 
lfci7/2 orbital in 209 Pb as an example. Although, on the level of one- and two-phonon admixtures, 
the fully chaotic GOE regime is not reached, the eigenstates of the model carry significant degree 
of complexity that can be quantified with the aid of correlational invariant entropy. With artifi- 
cially enhanced particle-core coupling, the system undergoes the doubling phase transition with the 
quasiparticle strength concentrated in two repelling peaks. This phase transition is clearly detected 
by correlational entropy. 

PACS numbers: 



I. INTRODUCTION 

It is now firmly established that many-body quantum 
dynamics acquire typical features of quantum chaos in 
the region of high level density. The interrelation between 
many-body physics and quantum chaos was studied, in 
particular, in the framework of the nuclear shell model 
with realistic and random interactions, see Q, 0, El> H 
and references therein. A description of a stable many- 
body system starts with the mean field and quasiparticles 
determined by the symmetry of the mean field. As exci- 
tation energy and level density increase, residual interac- 
tions convert stationary states of independent quasiparti- 
cles into exceedingly complex superpositions of many ba- 
sis states. The many-body structure at this stage should 
be described in statistical terms, and in the limiting case 
of extreme chaos it locally approaches the predictions 
of the Gaussian Orthogonal Ensemble (GOE) of random 
matrices. 

The complexity of stationary many-body states can be 
quantified with the aid of such characteristics as infor- 
mation entropy or inverse participation ratio |2,blU>ul 
of generic wave functions. These quantities smoothly 
change with excitation energy revealing strong mixing of 
original states and loss of simple quantum numbers, typ- 
ical features of quantum chaos. However such measures 
of complexity are unavoidably associated with a start- 
ing basis, in these cases determined by the mean field. 
They actually provide the localization length of a com- 
plex state in a chosen basis with no information on the 
correlations still present in the wave function. 

Different characteristics can be used in order to probe 
the sensitivity of the system to external perturbations. 
In classical case such sensitivity is known to be the main 
property of dynamical chaos. A special entropy-like 
quantity, called invariant correlational entropy (ICE), 
was suggested in Ref. 7] as a measure of complexity 
related to the response of a system to a random noise 



included in the Hamiltonian. This quantity is by con- 
struction invariant with respect to the basis transforma- 
tions. It also reflects the correlations and phase relation- 
ships between the components of the wave functions as 
far as they are revealed in the response of the system to 
a perturbation. Taking the strength of the interaction 
as a control parameter, one can clearly see the quantum 
phase transitions as it was demonstrated in the interact- 
ing boson model Q and in the evolution of pairing in the 
shell model @. 

The problem of the interaction of a quasiparticle with 
an even-even core is of considerable interest for nuclear 
physics. It was quite well studied for low-lying excited 
states where the interactions of the quasiparticle predom- 
inantly with quadrupole and octupole surface vibrations 
are important. At higher excitation energy, a large num- 
ber of other nuclear modes influence the damping of the 
quasiparticle motion. The mixing of the simple mode 
with the states of the next levels of complexity leads to 
the fragmentation of the single-particle strength over a 
wide domain of excitation energy, — the single-particle 
state obtains a spreading width. 

It was shown in Ref. [ljj that the spreading oc- 
curs in two stages. At the first stage of fragmentation, 
the single-particle state is spread over several doorway 
states. At the second step, the doorway states are spread 
through the mixing with many complex excitations re- 
lated to other degrees of freedom. The limiting case of 
self-consistent hierarchy of multi-step spreading was dis- 
cussed in Ref. [Tl| . 

The highly-excited single-particle mode has been ex- 
perimentally studied via one-nucleon transfer reactions 
[lOj. For example, broad structures located around 10 
MeV excitation energy have been observed in the re- 
action 208 Pb(a, 3 He) 209 Pb. It was shown ^ that the 
structures are connected with the excitation of high-spin 
orbitals IJ13/2, 2/i 11 / 2 and lfci7/2- The fragmented wave 
functions based on these states contain many compo- 
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nents. The related complexity of the wave function was 
analyzed in Ref. 0. It was concluded that the damp- 
ing of a simple mode cannot be unambiguously associated 
with the random contribution of more complex compo- 
nents to the structure of the states. 

The goal of the presented study is to introduce a new 
measure of the complex structure of highly-excited states. 
The paper is organized as follows. In Sec. II the for- 
malism of correlation entropy is introduced. In Sec. Ill 
the main features of the used model for the structure of 
excited states are presented. The complexity of of the 
excited states and the connection to the ICE is discussed 
in Sec. IV. The doubling phase transition in the regime 
of strong coupling in odd nuclei is demonstrated in Sec. 
V. Finally, in Sec. VI the conclusions are drawn. 



II. CORRELATIONAL ENTROPY 

To study the sensitivity of the excited states to varia- 
tion of external parameters we use the invariant correla- 
tional entropy (ICE) 0. The ICE method presumes that 
Hamiltonian H(X) of a system depends on a random pa- 
rameter A. The parameter A ("noise") is considered as a 
member of an ensemble characterized by the normalized 
distribution function V{X), 

J dXP(X) = 1. (1) 

In an arbitrary primary basis | k) , we follow the evolution 
of any stationary state |a; A) as a function of A. At a 
given value of A, the state can be decomposed as 

\a;\)=Y,CZ{\)\k), (2) 
k 

The ICE is defined as 

S a = -Tr{ Q a \n{e a )}, (3) 

where g a is the density matrix of the state |a) averaged 
over the noise ensemble. In the basis prior to the 
averaging we construct first this matrix for a given value 
of A as 

&,(A) = ^(A)Q*(A), (4) 
and then average over the ensemble, 

et k , = J d\P(\) e % k ,(\). (5) 

While the density matrix g a , eq. |@J, of a pure state, be- 
ing a projector onto the state \a) and having correspond- 
ingly only one nonzero eigenvalue equal to one, leads to 
zero entropy S 1 ", the averaged density matrix JSJ) has its 
eigenvalues between and 1 and produces non-zero cor- 
relational entropy. 



In contrast to information (Shannon) entropy I a of 
the same state |a), conventionally used for quantifying 
the complexity, 

7 « = _£|<7(«)|2 ln( | C («)| 2)) (6) 

k 

the ICE is basis-independent von Neumann entropy that 
reflects the correlations between the wave function com- 
ponents which are subject to fluctuations determined by 
the parameter A. The value S a for a given state typically 
increases with the complexity of the state and reaches the 
maximum at the point where the change of the param- 
eter around some average value implies the most radical 
change of the structure of the system. Such a point (in 
fact, a region in finite systems) can be identified with the 
quantum phase transition or crossover, and in the vicin- 
ity of this point the structural fluctuations of the wave 
functions are strongly enhanced. At further change of 
the average value of the noise parameter, the ICE, as a 
rule, goes down. This means that the stronger interaction 
has established a new order with greater rigidity with re- 
spect to fluctuations of parameters. In this way the sharp 
boundaries between the different symmetry classes very 
confirmed |8[ in the interacting boson models, and the 
critical strengths of isovector and isoscalar pairing in the 
sd shell model were established |g. 



III. THE MODEL OF EXCITED STATES IN 
ODD NUCLEI 

We adopted the Quasiparticle-Phonon Model (QPM) 
by Soloviev at al. [Tj| to describe the properties of highly- 
excited states in odd spherical nuclei. According to the 
model, the Hamiltonian of the system of an odd number 
A + 1 particles has the form 

H = h + H COIC + -ffcoupi- (7) 

The first term, h, describes the motion of a quasiparticlc 
in a mean field potential U created by the even-even core, 

h = -±Vi + U. (8) 
2m 

The core Hamiltonian is a sum of single-particle Hamil- 
tonians hi and two-body residual interactions Vij , 

A A 

H COIC =^hi + Y^ V id ■ ( 9 ) 

i=l i<j 

The last term in eq. Q is a sum of interactions between 
the odd quasiparticle and particles in the core, 

A 

ffcoupl = Yl V °>* ■ ( 10 ) 

t=l 
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In the spirit of the QPM, the Hamiltonian H corc is treated 
in the random phase approximation (RPA), i.e. the 
particle-hole configurations are built with the subsequent 
RPA diagonalization. 

The properties of the (A+l)-nucleus can be described 
in terms of the quasiparticle states a\ | 0) , quasiparticle- 
plus-phonon states [aj eg) Q\] | 0) and quasiparticle- 
plus-two-phonon states [a \ ® [Qj, ® QJ,]] | 0), where 
all combinations have the same total spin and parity 
quantum numbers J W M. Here a\ is the quasiparti- 
cle creation operator with shell-model quantum numbers 
a = (n,l,j,m), whereas Q\ = Qw denotes the phonon 
creation operator with the angular momentum A, projec- 
tion /i and RPA root number i. 

The following wave functions describe in the QPM the 
ground and excited states of the odd nucleus with angular 
momentum J and projection M: 

V}m = C?W JM + ]T D» a {Ja) [alQl] JM 
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FIG. 1: Distribution of the single-particle strength of 
the state lk 17/2 in 209 Pb. 



(11) 



The wave function \I>o represents the ground state of 
the neighboring even-even nucleus (also quasiparticle and 
phonon vacuum) and a stands for the number within 
a sequence of states of given J 71 ". The coefficients Cj, 
D^(Ja) and F" v (I; J a) are the quasiparticle, quasipar- 
ticle + phonon and quasiparticle + two phonons ampli- 
tudes, respectively, for the state a. The norm of the wave 
function reads: 

(*jmI*jm> 



= ra 2 i + E^( Jq )] 2 + 2 E [*rV; . 

y av aw' ) 

The Hamiltonian J7J contains several parameters. The 
mean field JSJ was chosen to be the Woods-Saxon po- 
tential. The residual interaction JJjJ in the particle- 
hole channel was taken in a separable form with in- 
teraction strengths in each mode considered as ad- 
justed parameters. The quasiparticle-phonon coupling 
term (|10() does not contain any additional freedom. 
The core excitations with all spins and natural parity, 
1 — , 2 + , 3~, 4 + , 5~, 6 + , 7~ and 8 + , were included in the 
calculations. For each momentum and parity, the RPA 
states up to 20 MeV excitation energy were taken into 
account. The large phonon basis is necessary for correct 
description of the single-particle strength distribution in 
a broad range of excitation energy. 

The used set of parameters have been successfully ap- 
plied to describe the properties of low-lying as well as 
highly-excited states in 209 Pb, see for example Ref. [l5| . 
For studying the ICE, we selected the single-particle 



1&17/2 state. The properties of this state are studied 
in detail in Ref. |l5|]. This orbital is quasi-bound in 
the Woods-Saxon potential being located at 4.88 MeV, 
energy much higher than the Fermi level of 209 Pb. Be- 
cause of its high energy, the state is surrounded by 
many quasiparticle-plus-phonon and quasiparticlc-plus- 
two-phonon states. 



IV. COMPLEXITY OF STATES AND 
CORRELATIONAL ENTROPY 

At high excitation energy, the level density is large and 
it is convenient to calculate the single-particle strength 
distribution by means of the strength function [T(i| us- 
ing an averaging Lorentzian function. The distribu- 
tion of the single-particle strength [coefficient C 2 of eq. 
Ijlip] of the lfci7/2-state is shown in Fig. 1. The num- 
ber of quasiparticle-plus-phonon components included 
in the wave function is 420, while the number of 
quasiparticle-plus-two-phonon components is 1116. The 
Lorentzian smoothing parameter was chosen to be 0.2 
MeV. 

It is seen from Fig. 1 that the single-particle strength 
is spread over a broad interval of excitation energy. The 
largest fraction, 81 % of the strength, is concentrated be- 
tween 5 and 13 MeV excitation energy, and the state ac- 
quires the large spreading width in this domain, I> = 1.5 
MeV (T* is calculated as the second moment of the distri- 
bution). It was shown |l0lll6| that the main mechanism 
leading to the spreading is the interaction of the particle 
with the vibrations of the even core. 

In the case of 209 Pb the lowest core excitation is 3^~ 
state. The component [lj'15/2 <8> 3j~]i7/2+ dominates in 
the structure of the lowest two excited states, where 8 % 
of the single-particle strength is concentrated. The com- 
ponent [2g 9 / 2 C3) 4^ ]i7/2+ also contributes to the structure 
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of these excited states. The contribution of the com- 
ponents, where the particle is coupled to 2f, 4^ and 
6^ phonons, is more important at the higher part of 
the distribution (above 13 MeV). The coupling of par- 
ticle with 5~ phonons is significant in the domain of the 
main peak [15l |. The 7~ and 8 + phonons mainly con- 
tribute at excitations around 10 MeV |15l| . More complex 
quasiparticle-plus-two-phonon components influence pre- 
dominantly the secondary fragmentation of quasiparticle- 
plus-one-phonon components. 

The correlation entropy connected with the excited 
states was calculated using two types of the distribution 
function, the Lorcntzian function, 



1 



27T A 2 + ^' 

and normal (Gaussian) distribution, 
V{\) = 



1 —A 2 
exp— . 



(TV27T 



2ct 2 



(13) 



(14) 



The value of the parameters used in the calculations was 
A = 0.8 and a — 0.5. For this choice of the parameters 
both distributions have the same maximum value. 

The calculation was done as follows. For the values 
of A = A c the matrix elements of 7J CO upi, eq. (|10|l . are 
multiplied by the corresponding value of V{X C )- The 
new matrix elements are used to calculate the eigenvalues 
and eigenvectors of the Hamiltonian Q). The elements 
of the density matrix g% k i(X c ) are calculated according 
to the eq. I@J and averaged, eq. ©. In the case of 
the Lorcntzian distribution, the value of A c is changed 
in the interval [—3 A, 3 A], while in the case of the nor- 
mal distribution function A c is taken from the interval 
[-2.5cr, 2. 5a-]. 

To estimate quantitatively the absolute value of the 
ICE we use also the summed correlational entropy of N 
states, 



N 



F(N) = J2 SC 



(15) 



where a labels the eigenstates in the sector with given 
quantum numbers. The normalized value of F(N) is de- 
fined as 



F 



1 

N 



F(N). 



(16) 



The wave function (|llfl includes many components. 
Because of the different structure of the components, the 
values of the matrix elements of H coup i reveal large fluc- 
tuations. There are many small matrix elements and the 
corresponding components are weakly correlated. Their 
contribution to the spreading process has to be negligible. 
The ICE gives the opportunity to construct the appropri- 
ate basis including only the most important components. 
The dependence of the ICE on the dimension of the basis 



Truncation [ % of max. matrix element] 

FIG. 2: Dependence of correlation entropy on the truncation 
of the quasiparticle-plus-phonon basis. 



is presented in Fig. 2. The basis is truncated according 
to the criterion connected with the maximum matrix el- 
ement. All components having the matrix element of the 
coupling with the original single-particle state less than a 
certain fraction of the maximum matrix element are not 
included in the calculation. Here only the one-phonon 
components [a\ ® Q\\ | 0) of the wave function l(TT|) are 
taken into account. It is seen that the normalized en- 
tropy, eq. I|16f) . reveals saturation for the small values of 
the truncation parameter. The components, whose ma- 
trix element is less than 5% of the maximum one, con- 
tribute weakly to F . The number of such components 
is large but their influence on the spreading process is 
negligible. 

The results for the ICE found with the two distribu- 
tion functions l|13|l and l|14f) are compared in Fig. 3. Fig. 
3a displays the correlation entropy of the obtained states 
in the case of the normal distribution, eq. I|14l) . There is 
a background of small values of the ICE due to the large 
amount of weakly interacting states. The greater val- 
ues of entropy are located in the vicinity of the peaks of 
the single-particle strength distribution pointing to the 
regions where the proximity of many strongly coupled 
states leads to enhanced sensitivity of the wave functions. 
As mentioned above, in the region around 8 MeV the par- 
ticle is coupled predominantly with the 3^ phonon, while 
around 13 MeV the coupling is occurring mainly with the 
first positive parity phonons. One can see the tendency 
to the enlargement of entropy when the excitation en- 
ergy increases and a larger density of the states at higher 
excitation energy enters the game. 

The calculation of the ICE with the Lorcntzian dis- 
tribution function shown in Fig. 36 leads to the results 
similar to those in Fig. 3a. The regions including the en- 
hanced density of states with large entropy appear at the 
same excitation energy. In the case of the normal distri- 
bution the value of F(N), eq. I|15|l . is 37.50, while for the 
Lorentzian it is 54.44, mainly because the Lorentzian has 
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FIG. 3: Correlation entropy of the excited states 17/2 in 
209 Pb. a) The normal (Gaussian) distribution function is used 
in the calculation for upper panel, and b) the Lorentzian dis- 
tribution is used for lower panel. The wave function includes 
only quasiparticle-plus-phonon components. 



the long tails. In spite of this difference in the value of 
F(N), the main features of both distributions of the ICE 
are similar being determined mostly by the interaction 
and density of states and only marginally by the type of 
the noise distribution function. 

To test the sensitivity of the results to the parameters 
of the model, we introduce a new parameter k. Multi- 
plying the matrix elements of i? CO upi, eq. I|10|l . by this 
number we vary the overall strength of the particle-core 
coupling. The value of k — 1 corresponds to the realistic 
strength of -ffcoupi chosen according to the initially fitted 
parameters. This case is presented in Fig. 3. Only the 
components [a\ (£> Ql] | 0) of the wave function are 
taken into account. 

The absolute value of entropy F(N), eq. Ijl5|l . depends 
on the value of k, i.e. F(N) F(N,k); the function 
F(N, k) is shown in Fig. 4 revealing a pronounced maxi- 
mum at k c = 1.6. The maximum of F{N 1 k) at k = k c can 
be identified with a quantum phase transition smeared 
as expected in a finite system. The single-particle state 



FIG. 4: Dependence of the integral correlation entropy on 
the strength of particle-core coupling. Only the components 
[c*l®Ql<] I 0) of the wave function lilt are taken into account. 



interacts strongly with the sea of more complex excita- 
tions, and its contribution to the structure of any indi- 
vidual excited state is strongly reduced. But it has to be 
pointed out that even for the case of maximum entropy, 
k = 1.6, the single-particle component preserves its dom- 
inance in the main peaks of the strength distribution. 
When the more complex quasiparticle-plus-two-phonon 
components are included in the wave function i|ll|l . the 
correlation entropy rapidly increases. The value of nor- 
malized entropy (eq. I16JI for the case presented in Fig. 
1 is larger than the maximum value shown in Fig. 2 
by a factor close to 6. The additional correlations are 
mainly due to the coupling of qusiparticle-plus-one- and 
quasiparticle-plus-two-phonon components. Because of 
this hierarchy of couplings, the distribution of the single- 
particle strength is not affected too much by the new 
terms in the wave function l|ll(l . 

For small values of k, the mixing of states is suppressed 
and the ICE decreases rapidly, see Fig. 4. This case 
corresponds to the weak particle-core coupling and rel- 
atively simple wave function. The strength of single- 
particle states is distributed in a narrow vicinity of its 
unperturbed energy. The nearest-neighbor level spacing 
distribution for k = 0.2 is very close to the Poisson distri- 
bution. In the case when only quasiparticle-plus-phonon 
components are included in the wave function (|llf> . even 
at k = 1 the correlations are quite large, Fig. 4. This 
influences the level spacing distribution of 17/2+ states. 
The nearest-neighbor spacing distribution of 17/2 + ex- 
citations can be fit by the Brody distribution |l| with 
the Brody parameter equal to 0.4. At this stage the 
single-particle component is not completely smeared and 
chaotically distributed over the excitations of the system. 
Clear remnants of the simple excitation seen in the com- 
plicated wave functions are similar to the phenomenon of 
scars |l7j| existing in simple quantum systems. 

The nearest-neighbor spacing distribution of 17/2+ cx- 
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FIG. 5: Nearest-neighbor spacing distribution of 17/2 states 
in 209 Pb. The calculated values are fit by the Brody distribu- 
tion with the Brody parameter equal to 0.6. 



FIG. 6: Distribution of the single-particle strength of the 
state 1&17/2 in 209 Pb for large value of particle-core coupling 
after the doubling phase transition. The value of k is k = 2.0. 



citations for the case when the quasiparticle-plus-two- 
phonons components are included in the wave function 
(|ll|l is shown in Fig. 5. The distribution is calculated for 
the region of the main single-particle strength, 5-13 MeV, 
where the density of the states is not changed much. 
As seen from Fig. 5, the distribution can be described 
by the Brody distribution with the parameter equal to 
0.6. The higher degree of chaoticity indicates the im- 
portance of more complex components and their influ- 
ence on the damping process. At the first stage, when 
the single-particle states interact only with quasiparticle- 
plus-one-phonon components, the regularities induced by 
the mean field, such as the structure of the level density, 
are not completely destroyed. At the next stage, when 
the interaction with the components of the next level 
of complexity is switched on, these regularities are partly 
smeared, but the wave function 1|1H is still relatively sim- 
ple. One can recall old results [23j, where the coupling 
of an unpaired particle with the collective monopole and 
quadrupole modes was considered in detail (later these 
results were applied [24[) to the spreading width of gi- 
ant resonances due to their mixing with low-lying shape 
vibrations). In exactly solvable models with a particle 
attached to a single level it was shown [2j| that the main 
effect of the particle-phonon coupling is in creating a co- 
herent state of the phonon field. The chaotic elements 
should be associated with the mixing between various 
quasiparticle levels and coupling with different phonon 
modes. 

More multi-phonon components are to be taken into 
account for reaching the full complexity. As degree of 
complexity grows, the wave function (|llf) approaches the 
realistic wave function necessary to describe the complex 
stricture of high-lying excited states. It has to be pointed 
out that a few quasiparticle components are weakly influ- 
enced by the growth of complexity. By this reason, even 
a rather simple wave function, truncated for example, by 



the quasiparticle and quasiparticle-plus-phonon compo- 
nents can be used to describe the transition probabilities 
connecting high-lying and low-lying excited states. For 
this purpose, a simple practical procedure of averaging 
high- lying single-particle strength can be used (2{|. 



V. DOUBLING PHASE TRANSITION 

For large values of k the particle-core coupling becomes 
very strong. The single-particle strength is split into two 
main pieces repelled to low and high excitation energies. 
The case for k — 2.0 is shown in shown in Fig. 6. The 
single-particle strength corresponding to the low and high 
energy peaks is 31 % and 21 %, respectively. The rest of 
the strength is distributed in between the peaks. 

To identify the physical nature of the quantum phase 
transition that occurs at an enhanced particle-phonon 
coupling, we can analyze the behavior of the strength 
function of the original single-particle state coupled with 
states of a more complex character. Such an analysis was 
first presented in Ref. |18| in application to giant reso- 
nances. Our formulation is in fact a particular limiting 
case of a generic problem [l9| of a "bright" state coupled 
to the background that contains a hierarchy of states of 
increasing complexity. In this consideration, the bright 
state |0), in our case the bare single-particle state with 
unperturbed energy Eq, is mixed with the background 
states that contain, apart from the quasiparticle, the 
phonon excitations of the core with unperturbed energies 
e„. The eigenstates \a) that emerge from this mixing are 
complicated superpositions 

| a )=C«|0)+]TC». (17) 
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Their energies E a are the roots of the secular equation 



VI. CONCLUSION 



X{E) = E-E Q ~J2 



IKI 



E-e„ 



0, (18) 



where V v are the matrix elements of coupling of the bright 
state with the background states. The fragmentation of 
the strength over the eigenstates \a) is given by 



dX 
~dE 



E=E„ 



E 



IK 



(E a - e u y 



(19) 



The strength function of the bright state is found as 



^H^TiCo"! 2 ^-^). 



(20) 



If the quantities V 2 , the average value (|K| 2 ) of cou- 
pling matrix elements squared, and D, the mean level 
spacing in the background, can be considered to be 
weakly fluctuating from one state \v) to another, the 
result of the mixing is determined by the ratio of the 
"standard" spreading width, 



r. 



2tt 



D 



(21) 



to the energy range a ~ ND, where N is a number 
of effectively interacting background states. Typically, 
one has in average a Breit-Wigner strength function 
with the width (|21(l for T s <C a; the further evolution as 
a function of increasing coupling leads to the Gaussian 
sha pe with the width increasing towards and beyond a 
pll2?il2lll22T|. In the transitional region the dependence 
of the strength function upon the coupling intensity V 
changes from the quadratic to linear |18L |20| . Finally, 
at even stronger coupling, eqs. (|18|) and l|19|l predict a 
phase transition to the new situation when the strength 
is accumulated in two peaks on both sides of the centroid, 
see Fig. 2 in Ref. 18]. In the extreme limit, the peaks 
arc pushed out of the region of the background states, 



E « En ± 



5> 



1/2 



(22) 



and the strength of the original state, according to eq. 
is evenly divided between the peaks. 
The physical mechanism of this doubling phase transi- 
tion (known also in quantum optics) is the following. In 
the regime of strong coupling, the background states of 
the same symmetry turn out to be intensely interacting 
among themselves through the bright state. Effectively 
this interaction is close to the factorized one built as a 
product of the matrix elements to and from the bright 
state. The factorized coupling leads to the formation of 
the second collective state as a coherent superposition of 
the background states. This state accumulates a signifi- 
cant strength, and the repulsion between the two collec- 
tive states leads to the two-peak pattern. We see that 
the doubling phase transition is adequately identified by 
the maximum of correlational entropy. 



In the presented study the recently suggested mea- 
sure of the complexity was tested for the wave function 
of a high-lying excited quasiparticle state in a system 
where the fermionic quasiparticle strongly interacts with 
bosonic collective excitations of the core. The new quan- 
tity — invariant correlational entropy — was used to es- 
timate the growth of complexity as a result of admixture 
of many new components to the wave function of a quasi- 
particle. 

In contrast to information entropy that displays the 
degree of complexity of individual states with respect to 
a certain reference basis, the ICE is basis-independent 
reflecting mostly the sensitivity of a given state to the ex- 
ternal noise. It was shown that the ICE depends mainly 
on the interaction strength and density of the background 
levels but less on the distribution function used as a noise 
generator. 

We have also calculated the neighboring level spac- 
ing distribution as a conventional indicator of quantum 
chaos. It is shown that this distribution is correlated with 
ICE moving in the direction of the Wigner distribution 
characteristic for the Gaussian orthogonal ensemble (but 
not reaching this limit) in parallel to the increase of the 
value of ICE. Although the model wave function trun- 
cated on the level of the quasiparticle-plus-two-phonon 
components is not sufficiently chaotic to manifest entirely 
the complicated structure of high-lying excited states, it 
can be used to describe the distribution of a few quasi- 
particle components. These components in fact play the 
role similar to the scars in quantum chaos. 

As a function of the overall strength of quasiparticlc- 
phonon interaction, the ICE increases and reaches a pro- 
nounced maximum at the coupling strength equal to 1.6 
of its realistic value. In this region the system would un- 
dergo the quantum doubling transition when the original 
single-particle excitation forms its own counterpart built 
of more complicated states with the same quantum num- 
bers and the two peaks repel each other and share the 
original strength. The ICE properly reflects this trans- 
formation. The question remains whether such a phe- 
nomenon, known in quantum optics, could be observed 
in real nuclear spectra. 

Similar calculations have been also performed for the 
excited neutron deep-hole orbital l/iii/2 in 207 Pb and 
high-lying proton orbital H13/2 in 145 Eu. The obtained 
results are in agreement with the above conclusions 
which confirms generic features of underlying physics. 
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